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We have obtained a semiclassical approximation of the matrix elements of a quantum propagator 
for coherent sates basis, for that propose we have used the symplectically invariant Weyl repre- 
sentation of the propagator and gone beyond the common stationary phase approximation. The 
comparison with a "realistic" chaotic system that derives from a quadratic Hamiltonian, the cat 
map, reveals that the expression here derived is exact up to quadratic Hamiltonian systems. 
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Path integrals appear as a useful calculation tool for 
many quantum and statistical mechanical problems 
while coherent sates are widely known to represent quan- 
tum states with the most classical resemblance, in the 
case of an harmonic oscillator they obey the classical 
equations of motion and are minimal uncertain states. 
Also, the coherent sates form an overcomplete basis 
which is a necessary ingredient for the construction of a 
path integral 0| . It implies the existence of several forms 
of path integrals, all quantum mechanically equivalent, 
but each leading to a different semiclassical limit. The 
formulation of path integrals applied to coherent sates 
has become widely used in many areas despite the very 
shaky mathematical background [3]. As was recently 
pointed out, the coherent states path integral breaks 
down in certain cases where neither the action sug- 
gested by Weyl ordering nor the one constructed by nor- 
mal ordering gives correct results when the Hamiltonian 
involves terms that are non linear in generators. In order 
to understand the quantum classical limit it is imperious 
to have a correct semiclassical expression of the quan- 
tum propagator in the most classical states, that is, in 
coherent states. In this Letter, by using a symplectically 
invariant formulation, we derive a semiclassical expres- 
sion for the coherent states propagator that goes beyond 
the common stationary phase approximation [5|. The 
comparison with a "realistic" chaotic system that derives 
from a quadratic Hamiltonian for which the semiclassical 
limit is exact, the cat map, reveals that the formulation 
is exact up to quadratic Hamiltonian systems. 

The coherent states \X) centered on points X — (P, Q) 
in phase space, are obtained by translating to X the 
ground state of the harmonic oscillator, its position rep- 
resentation is 



erality. The overlap of two coherent states is then 
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{X \X'} = exp 
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Where the wedge product 

X AX' ^ PQ' - QP' ^ (JX) .X', 

In what follows we use the centers and chords formulation 
developed by Ozorio de Almeida Q for both classical and 
quantum mechanics. The matrix elements of the unitary 
propagator J7* , that governs the quantum evolution of an 
L degrees of freedom system, in the coherent states basis, 
<Xi\U'\X2 >=U\Xi,X2) are 



U\Xi,X2)= i^—j J <X,\U\x)R,\X2>dx. (3) 

We have inserted the, symplectically invariant, Weyl de- 
composition of the propagator in terms of reflection op- 
erators [6] 
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dxU\x)R^, with U\x)=tr 



(4) 

with Rx the set of reflection operators thought points 
X — {p, q) in phase space 0| then the operator J7*is 
written in terms of the function in phase space, U^{x), 
its center or Weyl Wigner symbol. The action of the 
reflection operator R^ on a coherent state \X) is the x 
reflected coherent state ESl: 



RJX) 



XAx 



2x-X) . 
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For simplicity, unit frequency (w = 1) and mass (m — 1) 
are chosen for the harmonic oscillator without loss of gen- 
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So that, inserting ([5]) and ^ in Q the propagator in 
coherent states becomes, 



. 1 iXiAXo 

U\X^,X2) = -—j:e^^ 

{■Kh) 



(6) 

with ^0 = (-'^i — X2) the chord joining the coherent 
states points, and X = ^^^^'^ is the mid point. 
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For an Hamiltonian system, the propagator is = 
exp[— with H the quantum Hamiltonian operator 
of the system. As has been seen in @ , that for sufficiently 
small times the Weyl symbol of such a propagator can be 
written as 



U'{x) = |det[l-(a^-H)2]|3exp 



with H(x) — tr 



RxH 



(7) 

the center or Weyl symbol of the 



Hamiltonian operator H and where % is the correspond- 
ing Hessian matrix evaluated at the point x. In analogy 
to the classical theory, H{x) will be within 0{K) the clas- 
sical Hamiltonian Hc{x). Also, if H{x) is quadratic ([7]) 
represents a unitary operator . 

For longer times the composing an even number K 
of unitary operators for times e = t/ K is performed. In 
the Weyl representation such composition results [6] 



U\x) 



dx 



K 



K 
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Ak + i{x,xx,--- ,Xk) 



(8) 



with dxf = dxK'-'dxi and Ai<-+i(x, xi, • • • ,xk) the 
symplectic area of the polygon with endpoints centered 
on x and whose j'th side is centered on Xj. This formula 
is exact, while in the limit K ^ oo, the expression ([7]) 
can be inserted for Ut_ {x). In which limit the amplitude 
vanishes, yielding to the path integral 



U\x) = lim 



dx 



K 
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KL 
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Analogously, for the case of classical composing of an 
even number K flows for time e = t/ K in the limit where 
K ^ oo^ the center generating function of the canonical 
transformation is 

t ^ 

S\x) = lim -—y^Hc{xj) + AK+iix,xi,---,XK)- 

A'— !-OC K ^ — ' 

(10) 

The polygon Ak+\ has one large side ^ passing through 
the center x and K small chords tangents to the orbit 
as X oo such that dS^ jdxj = 0. Hence, the center 
variational principle is obtained: The center action 



S'^ix) = (j)^p ■ dq - J H,{x)dt 



(11) 



is stationary along the classical trajectory 7. The paths 
to be compared always have their endpoints centered on 
the point x. The second integral is evaluated along this 
path, whereas the first integral is closed off by the chord 
centered on x. 

Hence, the phase of the integral in ^ coincides with 
the center action S''(x) for the polygonal path 7. The 
center variational principle ensures that this center ac- 
tion is stationary for the classical trajectories centered on 



X. The semi-classical approximation for the Weyl prop- 
agator was obtained by performing a stationary phase 
approximation in Q then the center variational princi- 
ple yields @ : 



U\x) 



sc 



E 



2^cyi^{ih-^Si^{x)+ia\] 
|det(X.., + l)|^ ' 



(12) 



the summation is over all the classical orbits 7 whose cen- 
ter lies on the point x after having evolved a time t Q. 
Where S** (a;) is the classical center generating function of 
the orbit, A^^ is the symplectic matrix for the linearized 
transformation between the neighborhood of the tips of 
the chord ^-y joining the initial and final point of the or- 
bit. For sufficiently short times such that the variational 
problem has an unique solution there will have a single 
chord. Although for longer times, there will be bifurca- 
tions producing more chords. In the case of a single orbit, 
the corresponding Morse index al^ — 0. 

In order to obtain a coherent state path integral, ex- 
pression ^ is inserted in the coherent state propagator 
del) so that 



U\Xi,X2)^ lim 



K- 



dx 



KL 
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+\k+i{x,xx,--- ,xk)\- 



Recalling the linear relation of A^+i with x , 



(13) 



K 

Ak+i{x,xi, ■ ■ ■ ,xk) = C + ^ a X, with^ = ^(-l)-'^^. 

(14) 

The X integral in (|13|) is quadratic and can then be solved 
exactly yielding 



. iXiAXo 

U\Xi,X2) = e — lim 

K^o 



dxf 1 t \2 



xe 



i{AK+i(^,a;i,---a;/f )-i E^=i ^(^")}(15) 



where ^ is the chord passing through the mid point X. 
Note that according to the chord ^ depends only on 
the other centers xi,--- ,xk ■ The usual semiclassical 
limit of the coherent states propagator is obtain, in the 
limit where fi — > , for a stationary phase approximation 
of ([15]). The phase in (|15p is the center generating func- 
tion pop . Hence, according to the center variational prin- 
ciple the path integral is stationary for classical orbits, 
then analogously to (|12p stationary phase approximation 
yields 



U\Xi,X2)sci^Yl 



2^e 



|det(X^ + l)|2 



(16) 
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where the index runs over ah the contributing classical 
orbits 7 with center generating functions'^ (X) whose cen- 
ter lies on the point X and is the chord joining the 
initial and final point of the orbit. Again, will have a 
single chord for sufficiently short times. However, as a 
consequence of a Gaussian cutoff on the length of the 
chords, the amplitude will be severely damped if the 
classical chord is long. We have define the action 
S^{x) = S!^{x) + hai^ in order to include the Morse in- 
dex in the action. Equation (|16p has been obtained by 
an equivalent procedure by dos Santos and de Aguiar [SJ . 
Also, as will be shown in Fig 1. expression ([T5)) does not 
give accurate results for quadratic Hamiltonian systems. 

An other way to obtain a semiclassical approximation 
for the coherent sates propagator is to insert the semi- 
classical approximation for the Weyl propagator in 
® so that 



^ , !^Pzil^z4^eM^^(-)+«oHd. (17) 
|det {M-y + l)\^ 

This procedure is equivalent to perform in (|13p the sta- 
tionary phase approximation in the dxn ■ • ■ dxi variables 
before integrating in dx. We must now perform the x in- 
tegral, in a first place let us perform an ordinary station- 
ary phase approximation in pzp . Expanding the phase 
of (jl7p around the stationary point a;o a-i^d performing an 
usual stationary phase approximation, another semiclas- 
sical approximation for the propagator in coherent states 
rises 



starts near X2 and ends near Xi will have an important 
contribution, this orbits will have a center point close to 
X. Hence, in order to perform the phase space integral 
in (jl7p . let us expand the center action up to quadratic 
terms near the mid point X so that 

Si^{x) = Si^{X)+l^x' + x'Bt^x' + 0{x'^) (19) 

with X = X + x' and 5'*(A) is the action of the orbit 
through the mid point X , while the symmetric matrix 
Bt^ is the Cayley representation of M.^ 



1 + M-, 



I d^S'^ix) 
2^ 



dx"^ 



(20) 



After the linearization of the flux around the mid point 
(fT9)) is used in ([T7)) a quadratic integral is obtained that 
can be exactly solved. Using 



V(detT^ 



(21) 



the semiclassical propagator in coherent states, takes the 
form 



\dct[Vt(M-, + l)]\2 



X exp ■ 



-\Xi AX2 + SUX) + 5Bt5 + M 



(22) 



with the complex matrix Vt — 1—iBt-y, such that det Vt = 
I det Vt|e'^, the point shift S = y2(<^— ^0), while the real 
matrices Ct and Bt are defined such that 



Ct~iBt = {l + iJBt^jy 



(23) 



U\Xi, X2)sc2 — 



^{i[s;«o) + iXiAX2]-i(xo-X)'} 



|det(X^ - 1) 



where a;o must satisfy ^(xq) 



J 



dx 



(18) 



that is the chord through xq equals the chord joining the 
coherent sates, and the chord action ^'(^o) — S^i^o) — 

^0 A xq [6] being the Legendre transform of S"* (xq). This 
expression (fT8)) is complementary to equation ((T6)) . it is 
written in terms of chords instead of centers. Although, 
as is the case for the chord action the expression 
(lisp diverges for very short times where the monodromy 
matrix becomes the identity. 

In obtaining (jlSp it was assumed that the Gaussian 
term in (|17p is smooth close to the stationary point. Al- 
though, there is a compromise between the point of sta- 
tionary phase and the maximum of the Gaussian in X 
. If a stationary phase point falls far from X its con- 
tribution will be severely damped. For this reason, it is 
not surprising that expression psp fails, as is shown in 
Fig 1. Also, it must be noted that classical orbits that 



The expression ([22t of the semiclassical matrix ele- 
ments between two coherent states of the quantum prop- 
agator, is the main contribution of this Letter. It is ex- 
pressed entirely in terms of classical objects, namely the 
action S** (X) of the classical orbit whose mid point is X 
, the point shift S, the monodromy matrix A4j and its 
Cayley representation Bfy . It is important to remark the 
Gaussian term that dampens the amplitude for large val- 
ues of the point shift S, that is for orbits that starts far 
from the points X2 (or end far from Xi). So as, the main 
contribution in the sum over classical orbits in will 
come from the orbits 7 , whose initial point Xi lies the 
closest to X2- It is important to remark that, the same 
expression ([^^ could have obtain by perform in ([TS]) the 
complete quadratic integral instead of a stationary phase 
approximation that led to (jl6p . 

Also, note that, if t = 0, the quantum propagator is 
just the identity operator in Hilbert space, the classical 
symplectic matrix is the identity, the center action is null, 
and the result ^ for the overlap of coherent sates is 
recovered 

We will study the matrix elements on a surface of 
section that is transversal to the flux, in analogy with 
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classical Poincare surfaces of section. The study of au- 
tonomous fluxes through a map on surface of section is 
a standard procedure, and quantum surface of section 
methods are shown to be exact fld\ for general Hamilto- 
nian systems. Here, the theory is applied to the cat map 
i.e. the linear automorphism on the 2-torus generated by 
the 2x2 symplectic matrix M, that takes a point x- to 
a point a;_|_ : x+ = A4x- mod(l). With the symplectic 
matrix 



M 



2 3 
1 2 



(24) 



This map is known to be chaotic, (ergodic and mixing) 
as all its periodic orbits are hyperbolic, it corresponds 
to viewing stroboscopically the motion generated by a 
quadratic Hamiltonian [7]. Quantum mechanics on the 
torus, implies a finite Hilbert space of dimension — 
, and that positions and momenta are defined to have 
discrete values in a lattice of separation 0, [ll| . The 
cat map was originally quantized by Hannay and Berry 
[ll| in the coordinate representation the propagator is: 



the exact expressions obtained with the quantum prop- 
agator on coherent states (f^^ . It can be seen that 
neither (ITSl) nor (fT5|) are good approximations of the ex- 
act coherent sates matrix elements giving errors in the 
amplitude of more than 10% or that highly grows with N 



□ 




< qfe|UAi|qi > = 



e « 



(25) 



where the states < q\qj > are periodic combs of Dirac 
delta distributions at positions q = 3/Nmod{l), with j 
integer in [0, — 1]. Bold fonts characters are used to 
describe quantum objects on the torus. The definition of 
the periodic coherent state [9], in accordance to ^ is 



< X|qfc > = 



j = -oo 



=.-M*-P(j+T-V«)+^(j+Q-V«)" 



(26) 



Indeed the construction on the torus from the plane is 
obtain in terms of averages of translations [0| . 

Figure 1 shows the relative error E of the amplitude 
of the semiclassical approximations U*{Xi, X2)sci Stiid 
U*(Xi,X2)sc2 obtained in ([TO)) and ^TE\i with respect to 



Figure 1: Relative error of the amplitude of the Semiclassical 
expressions U^{Xi, X2)sci and U^{Xi, X2)sc2 as a function 
of A'^. In dotted line error of U*{X\,X2)sci and in full line 
the error of (7'(Xi, ^2)502. 



To conclude, in this work the semiclassical coherent 
state propagator has been obtained using the Weyl repre- 
sentation. The expression of U*{Xi, X2)sc3 obtained in 



goes beyond a stationary phase approximation and 
is exact for linear systems. Also, the comparison with 
a system whose semiclassical limit is exact has allows to 
correctly check the exactness of the obtained expression 
up to quadratic Hamiltonian systems. 
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